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Abstract 

We construct and study a discrete time model describing the conflict in- 
teraction between two complex systems with non-trivial internal structures. 
The external conflict interaction is based on the model of alternative inter- 
action between a pair of non- annihilating opponents. The internal conflict 
dynamics is similar to the one of a predator-prey model. We show that 
the typical trajectory of the complex system converges to an asymptotic 
attractive cycle. We propose an interpretation of our model in terms of 
migration processes. 
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1 Introduction 

Since the beginning of 20-th century the Lotka-Volterra model of prey-predator 
interaction is one of the main models for simulation of many processes in popu- 
lation theory and economics. As a rule, continuous models where Lotka-Volterra 
equations have ratio-depended parameters are studied (see, for example [3J El El 
ITTj [121 [131 [T6| HU [SO]). Logistical and Ricker's models are also studied in some 
works, for example [5]. In the majority of works the prey-predator interaction is 
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treated only inside of single region, and no migration from one region to another 
is considered. 

In some works 0, [5] models with migration are studied, with a process of 
migration. 

There are also only few works (see [5] and references wherein), in which dis- 
crete models are considered, though in reality such processes are more natural, 
since they take better into account seasonal phenomena (reproduction, migration, 
etc.). 

The main aim of the majority of works is determination of stable points, 
bifurcation points, asymptotic behavior, and analysis of model's depending on 
the coefficients of the equations. 

In [3] synchronization of population dynamics with natural phenomena (like 
change of seasons and floods) is studied. In the work [20] dependence of pop- 
ulation dynamics on population density and migration is studied. In the work 
[6] migration is not assumed to be random, but aims at maximization of some 
function of the population. At last, in [3] the influence of stochastic terms in a 
Lotka-Volterra model is described, and interesting figures are presented. 

In recent works [T71 Salam and Takahashi study conflict models, similar to 
ours. They introduce important and more complex multi-opponent systems. In 
[TT] not only conflict, but cooperation between opponents is studied. The figures, 
obtained by them, are very similar to Figure 3 of the present work. 

In this work we construct a model that joins two most rarely studied variants 
of Lotka-Volterra model, i.e., a discrete model with migration. Here individuals 
migrate not randomly, but according to strategies, discussed in section 5. 

We construct the model of the conflict interaction between a pair of complex 
systems A and B. The system is a finite set of positive numbers: P = (Pi, . . . , P/v) 
for A and R = (Pi, . . . ,Rn) for B, where N means the quantity of parame- 
ters that characterize the system. We study dynamics in the discrete time. So, 
the evolution of every system is described by the sequence of vectors with non- 
negative coordinates P n = (pf n) , . . . , P$°) for A, and R n = (R { ™\ . . . , R^ ] ) for 
B, n = 1, 2, . . .. The vectors P and R correspond to the moment n = 0. Natu- 
rally, each system tries to reach the optimal values of its coordinates. In reality, 
due to the conflict interaction, every coordinate changes in a complicated way. 
The evolution of all changes is determined by double dependence: by the con- 
flict interaction between systems (which we shortly describe in section 3), and by 
the mutual "fight" of coordinates (of the prey-predator type interaction) inside 
every system. We suppose that every system is complex in the sense that its 
elements may be treated as one of the types: dominant (predators, employers) 
or dependent (preys, workers). So, every coordinate P/ n \ R^ may be regarded 
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as the quantity (population) of dominant, respectively dependent species at the 
position % at time n. 

The law of evolution inside of each (independent) system is described in section 
2. We suppose this law is identical in every system and is based on the well-known 
discrete Lotka-Volterra equation. 

In section 3 we shortly call the main results on conflict interaction between 
non-annihilating opponents. 

In section 4, that includes the main results of the work, we construct a dynam- 
ical system describing simultaneous conflict interaction both inside every system 
and between the systems. The outer interaction is an alternative conflict be- 
tween non-annihilating systems, whereas the inner one is a prey-predator model 
of Lotka-Volterra type. 

We may join these two types of interactions in a discrete time. Thus, our 
dynamical system consists of a discrete sequence of states. Two operations hap- 
pen at any fixed moment of time: redistribution of probabilities to occupy some 
controversial positions by opponent systems, and quantitative changes (namely 
population) of all species inside both systems. 

The computer modelling of such a complex interaction shows some very in- 
teresting phenomena. In this work we limit however ourselves to present only 
one observation. Namely, under an appropriate choice of parameters and initial 
data the complex system oscillates. We find a rather wide range for initial data 
for which the population trajectory in phase-space becomes cyclic. Moreover, we 
observe the stability of the limit cycle, so it is an attractor. 

2 Traditional models of population dynamics 

Malthus proposed in 1798 the population equation 



where P is the cumulative number of individuals (species), and b, d are the natural 
birth and death rates. In reality, one expects exponential solution 



describes the ideal population of biological species. The exponential rise, if b—d > 
0, or decrease, if b — d < at most in a local period of time. 

Verhulst introduced in 1838 a more realistic equation with saturation terms: 



dP 
~dt 



(b - d)P, 



(2.1) 



P(t) = P(0)e (b - d)t 



dP_ 
~dt 



{b - d)P - cP 2 , 



(2.2) 
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Figure 1: A typical shape of the logistical curve ^ = {b — d)P — cP 2 . 

where the coefficient c > represents the competition activity of individuals for 
living resources. The square power corresponds a conception of an alternative 
law of access to the living resource. 

The solution of (2.2) describes the 5-shaped logistical curve (see Figure 1) 
and corresponds better to the actual behavior of many population processes. 

The curve starts with a small value -P(O), exponentially increases, and then 
saturate at the capacity P(s) = — , b — d > 0. 

In the economic context, equation (2.2) can be written as follows 

^L = ( g -l)M- fM\ (2.3) 

where M is the capital (money), g and / are the average gain and loss percentages 
on the capital, and / stands for the coefficient of confrontation between individ- 
uals. If g — I < 0, the capital decays to exponentially; if g — I > 0, at the 
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beginning the capital increases exponentially quick, but then the growth slows 
down, so that it never reaches asymptotic value of saturation M(s) = s y L . 

Lotka (1907) and Volterra (1901) extended the Verhulst logistical equation to 
the Lotka- Volterra equations intended for the description of amount changes in 
populations of two species in interaction. These equations are also known under 
the name of predator-prey model. We will refer on Lotka- Volterra equations in 
the following form: 

^£ = aP- bPR - cP 2 
f = -dR + ePR-fR 2 , [2A) 

where all coefficients are nonnegative. 

The population of prey is described by the first equation. Without presence 
of any predators it grows exponentially at the beginning and then comes to the 
fixed capacity P(s) = a/c. The predators, without any prey to feed on, die out. 
When both species are present, the growth of the prey is limited by the predators, 
due to the term —bPR, and the predators grow if the amount of prey available, 
i.e. if ePR is large enough. 

There are many publications devoted to the analysis of Lotka- Volterra equa- 
tions (2.4) (see for example (16] and references wherein). 

The models with discrete time are also studied. In this case, equations (2.4) 
have the following view: 



p [n) = p in-l) + p(n-l) (a _ fep (n-l) _ c p(n-lh 
pin) = pin-l) + p 2 (n-l)(_ d + e p(»-l) _ fp^ 



(2.5) 
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Typical behaviour of discrete Lotka- Volterra model is shown in Figure 2. 



3 Conflict interaction between non-annihilating 
opponents 

In this section we shortly remind an alternative approach to describe the redistri- 
bution of conflicting positions between two opponents, say A and B, concerning 
an area of common interests. 

We consider the simplest case where the existence space of common interests 
is a finite set of positions Q = {ux, . . . , uj^}, N > 2. Each of the opponents A and 
B tries to occupy a position uii, i — 1, . . . , N with a probability P^(^j) = Pi > 
or Pb{uJi) = Ti > 0. The starting distributions of A and B along Q are arbitrary 
and normed: Ylf=iPi = 1 = Ylf=i r i- A and B can not be present simultaneously 
in a same position cjj. The interaction between A and B is considered in discrete 
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P[n]=P[n-1]*(a-b*R[n-1]-c*P[n-1])/500 
R[n]=R[n-1]*(-d+e*P[n-1]-f*R[n-1])/500 




Figure 2: Lotka- Volterra model with discrete time 



P[ n) =Pl n -'\a-bP{ 



(n-l) 



,(n-l) 



(n-l)> 



P 



(n-l), 



-d + eP 



(n-l) 



(n-lh 



a = 0.2, b = 0.006, c = 0.002, d = 0.008, e = 0.002, / = 0, 



PI 



(0) 
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time t G No- We introduce the noncommutative conflict composition between 
stochastic vectors p° = (pi, . . . ,px),r° = (r%, . . . ,r N ) e R+: 

p 1 := p° * r°, r 1 = r° * p°, p° = p, r° = r, 

where the coordinates of p 1 , r 1 are defined as follows 

<„ pf'd-arf) a, = rf>(l-«rf°>) 

where the coefficient -1 < a < 1, a / stands for the activity interaction. At 
the nth step of the conflict dynamics we get two vectors 

p n = p™" 1 * T n ~ l = p° * n r °, r n = r n_1 * p n_1 = r° * n p° 
with coordinates 

(n— l)/ n (n— 1)\ (n— 1)/-, (n— 1)\ 

(n) gj y (l-ar t v ; ) (n) rj ; (l-ap t v ; ) 

with z n a normalization coefficient given by 

z n = l-a{ V n -\v n - 1 ) 1 

with (•, •) the inner product in M. N . 

The behavior of the state {p n , r n } at time t = n for n — > oo has been investi- 
gated in [TJ IU d El El [10]. We shortly describe the results. 

Theorem. 1. For any pair of non- orthogonal stochastic vectors p, r 6 R^, 
(p, r) > ; and /ixed interaction intensivity parameter a ^ 0, — 1 < a < 1, 
wii/i condition a ^ J^Trj> ^ e sequence of states {p n ,r n } tends to £/ie /zrmt state 
{p 00 ,!- 00 } 

p°° = lim p n , r°° = lim r n . 

n—>oo n— >oo 

TTws /imii stoto zs invariant with respect to the conflict interaction: 

p oo = p oo^ r oo r oo = r oo^ p oo 

Moreover, 

' p°° 1 r°°, i/ p ^ r and < a < 1 

We emphasize that in the case of a purely repulsive interaction, < a < 1, if 
the starting distributions are different, then the limiting vectors are orthogonal. 
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3500 4000 4500 5000 



500 1000 1500 2000 2500 3000 3500 4000 4500 5000 

Figure 3: a = 1, p° = (0.5; 0.3; 0.2), p° = (0.48; 0.34; 0.18) 
p°° = (0.33; 0; 0.67), r°° = (0; 1; 0). 
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Therefore each of the vectors p 00 ,! -00 contains by necessity some amount of zero 
coordinates on different positions cu,. For example the typical limiting picture for 
p n , r n e M.+ is presented in Figure 3 (comp. with [T7J [T9]). 

If we start with a pair of identical vectors, p = r, then p°° = r°° too. That 
is, all non-zero coordinates of the limiting vectors are equal. 

In the general case, p, r 6 the coordinates p\ , r\ have at most sev- 
eral oscillations and then reach monotonically their positive or zero limits. The 
limiting values p^ ,^ 00 may be described in terms of starting states. 

Given a couple of stochastic vectors p, r e E" , p^r, (p, r) > 0, define 

D + := di, D- : = ^ di, 

where 

di=Pi- n, N+ := {i:di> 0}, N_ := {i : d % < 0}. 

Obviously 

< D + = -£>_ < 1, 
since p^r, and J^p* - £^ = = D + + D_. 

Theorem 2. Lei p/r, (p, r) > 0. In i/ie purely repulsive case, a = 1, the 
coordinates of the limiting vectors p°°, r°° have the following explicit distributions: 

f di/D, i6N + r °° = J ~ d il D i i£ N - (3 2 ) 

0, otherwise ' * ^ 0, otherwise, 

where D := D + = —D-. 

Remark. From (13. 2p it follows that any transformation p, r — > p', r', which 
does not change the values di and D, preserves the same limiting distribution as 
for the vectors p 00 ^ 00 . A class of such transformations may be presented by a 
shift transformation of coordinates, Pi — > p' { — Pi + a^, — > r[ = r, + aj with 
appropriated a^s. 

In the case — 1 < a < of the pure attractive interaction we have another 
limiting distribution. 

Define the set § := = r fc° = 0} and set 

S°°:={l,...,iV}\So. 

Theorem 3. In the purely attractive case, a = —1, the limiting vectors 
poo^oo are e g Ua i an( i fheif coordinates have the following distributions: 

p°° = r °° = J 1 / m ' % G § °° (3.3) 
2 1 0, otherwise, 
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where m = |S°°| denotes the cardinality of the set 

In general, it is an open question to give a complete characterization of So- 
Below we present several sufficient conditions for k to belong to the set §o- 

Simultaneously these conditions give some characterization for the points to be 

in 

We will use the following notations: 

a . ■= p . + n, Pi := p^i, <t\ := p\ + r\ p] := p\r\. (3.4) 



Proposition 1. If 

°i > o~ k , Pi > Pk, or o { > a k , p { > p k , (3.5) 

then 



PT = rT = 0, 



and therefore k G So- 

Proof. By (13.41) we have 



°l = Pk + r l ~ l/ z (pk + r k + 2p k r k ) = l/z(<r k + 2p fc ) 

where we recall that z — 1 + (p, r). Therefore each of the conditions (13. 5p implies 
that a} > a\. Further, since 

pi = l/z 2 (p k + (p k ) 2 + pkVk), (3.6) 

again from (I3.5p it also follows that p\ > p\. Thus, by induction, of > and 
pf > $ for all N > 1. 
Or, in other words, 

1 N 
1 Pi Pi Pi 

1< — < ^ < ... < -k-.., 

Pk p k Pk 

1 N 

K^<^<...<^...,iV = l,2,... (3.7) 

Pk n r% 

Thus, sequences of the ratios 

pf_ r[_ 
Pk' r k 
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are monotone increasing as iV — > oo. Assume for a moment that there exists a 
finite limit, 

1 < hm % = — = — = M < oo. 

N->°op% p™ pf l + r~ 

This is only possible if r°° = r£°, which contradicts (13. 7p . Thus, M = oo and 
therefore = 0, as well as r£° = 0. □ 

Let us consider now the critical situation, when for a fixed pair of indices, say 
i and k, the values Ok — Oi, pk — Pi have opposite signs, for example, Ok — &i > 0, 
Pk — Pi < 0. In such a case it is not clear what behavior the coordinates pf,rf 
and i r k w * n have when N — > oo. We will show that the limits depend on 
which of the two values, 2pi + Oi or 2p k + Ok, is larger. Moreover we will show 
that even if p^ is the largest coordinate, it may happen that = 0. Let for 
example, pk = maxj{pj, rj} and Ok = Pk + > Pi + 1 = however the value 
of Tk is such that = PkTk < Pi^i = p%- Then under some additional condition 
it is possible to have = 0. In fact we have: 

Proposition 2. Let for the coordinates Pi, r^p^, r^, i ^ k, the following 
conditions be fulfilled: 

o k > Oi (3.8) 

but 

Pk < Pi- (3.9) 

Assume 

2pk + o k <2p i + o i . (3.10) 

Then 

PT = rT = 0, (3.11) 

i.e., k G So 

Proof. We will show that (j3~^j) . f[3T9l) . and (ETTOj) imply 

p l k + r l k = ol<ol=p\+rl (3.12) 

and 

= Pk < Pi = Pi r l (3-13) 

Then (13.111) follows from Proposition 1. In reality (I3.12p follows from (I3.10P 
directly, without condition (13.91) . So, we have only to prove (I3.13p . 

With this aim we find the representation for p\ in terms cr, and o\. Since 
o\ = \jz(oi + 2pi) we have 

Pi = l/2{zo]-o i ). (3.14) 
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By (|£BH and (IXTIj) we get 

Pi = l / Z \p* + Pi + Pi°i) = - + 1 /2(^- " <Ti) + *i] 

= ^ Z °\ ~ ^)( 2 + Z °\ + *i) = J^ Z °\ + + Za \°* - 2(T * - ZC7 i a i ~ a l\ 

Therefore 

pi - p\ = 1/Z 2 [p k (l + p k + (T fc ) - Pi{l+Pi + Of)]. 

Thus, we have 

- p! = - ^) + ^) 2 - (^) 2 ) + ((^) 2 - K) 2 ) + 2fa - a*)] < o 

due to starting condition (13.91) . and (13.121) . Thus p\ < p], i.e., (13. 13[) is true. □ 

We stress that ( 13. lip is true in spite of > 0{. Of course, if au < cr« and 
Pk < Pi, then ( 13.111) holds without any additional condition of the form ( 13.101) . 

4 Model of conflict interaction between complex 
systems 

In this section we construct a dynamical model of conflict interaction between a 
pair of complex systems. Each of the systems is subjected to the inner conflict 
between their elements. For simplicity, we assume both systems to be similar 
and described by discrete prey-predator models of type (2.5). We introduce the 
conflict interaction between these systems using an approach developed in [fl [2j [U 
13 El El [10] . With such a rather complex situation we may obtain a wide spectrum 
of evolutions. In this work we study qualitative characteristics of the behavior of 
corresponding dynamical systems for some choice of parameters a, b, c, d, e, /, a 
(see (2.5), (3.1)) and values of initial populations of species Pi,Ri. 

The coefficient a, that shows intensity of the interaction between systems, has 
an important effect. The increasing a from zero to unit causes the appearance 
of a series of bifurcations. For a = we have two copies of independent Lotka- 
Volterra models. For small values of a both systems behave like pure Lotka- 
Volterra systems, coming them to a stable state. 

Under fixed parameters and the starting coordinates a = 0.2, b = 0.006, c = 
0.002, d = 0.008, e = 0.002, / = 0, P$ 0) = 3, P 2 (0) = 10, Rf ] = 5, = 20 we have 
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first bifurcation point at a ~ 0.0056781739. The coordinates (R\ ) oscillate 
and a cycle of a small period appears. 

The following increase of a shows the appearance of new bifurcation points 
that are characterized by an increasing value of the cyclic period. For the value 
a = 0.4815545975 a cycle of infinite period appears. This means that all coor- 
dinates rapidly reach the stable state. In this case some species may disappear, 
even if they had some stable positive values in a pure (a = 0) Lotka-Volterra 
model. 

The role of the coefficients a, b, c, d, e, / and initial quantity of the species 
Pi,Ri in a pure Lotka-Volterra model is well-known and described (see, e.g., 
[T6"l [12]). Partially, coefficients a,d govern the increase of the pray population 
when predators are absent and the predator population decreasing when prays 
are absent. In turn, the coefficients b, e are responsible respectively for the pray 
quantity decreasing with an increasing number of predators, and increase of the 
predator population with an increase of the number of prays. The last coefficients 
in each of the equations give the limitation of increasing of both populations. In 
other words, each population "makes pressure" on itself, it does not permit an 
infinite reproducibility. 

Questions about stable points, orbits, asymptotic behavior of orbits are well 
described for the classical Lotka-Volterra model. We shall recall that usually 
there are at least three equilibrium points. They are refereed in literature as 
follows (see, e.g. [16]): 



An equilibrium point is called stable point if after a sudden change of popu- 
lation it comes back to an equilibrium point some time later. This may happen 
monotonically, or with some oscillations. 

We should note that under the existence of stable points the behavior of the 
system is well defined by coefficients a, b, c, d, e, /. But under the absence of stable 
points, the behavior of the system is defined by the initial data Pi, Ri- Depending 
on how close the initial data are situated with respect to the equilibrium point, 
the system may evaluate in a different way. 

The role of all these coefficients is preserved in the case of our model. But 
now their influence is much more complex. We present here only first steps in 
this direction. We shall discuss not only stability zones, as it was pointed above, 



(1) trivial (0,0); 

(2) axial (a/6,0); 

(3) inner positive 





13 



but also the values of the coefficients for which the system oscillates along some 
closed cycles. 

The state of our dynamical system is fixed by a pair of vectors P n = (P[ n \ . . . , PjJ ), 
TL n = (bP,...,B%)) with non-negative coefficients, where n = 0, 1, . . . denotes 
the discrete time, N > 2 stands for the number of conflict positions. Here we 
study the most simple situation, when every system consists of only two agents: 
pray and predator, i.e. N = 2. The complex conflict transformation is denoted 
by the mapping 



where F is the composition of four operations, the specific mathematical trans- 
formations: F = \J\f~ 1 * Af]U. 

Let us describe them in an explicit form for the first step. 

The first operation U describes the interaction between elements inside every 
system separately according to the pray-predator model. Corresponding math- 
ematical transformation of vectors (the interaction composition) {P°, R } — -> 
{P°, R } is described by the system of equations of the form (2.5): 



where the passage to new values of coordinates is pointed by tilde, but not by 
changing of upper index, likely to (2.5). 

The following operation involves the interaction * (see (3.1)) between previous 
systems according to the theory of the alternative conflict for non-annihilating 
opponents (see, e.g. [U El HI CO El El EE])- To describe this operation we at first 
have to normalize the vectors P° = (Pi°\ P2 )i R° = {R\ ■> ^2 )j i- e -' to work 
with stochastic vectors. 

We use the following notation for normalization: A/"{P°,R } = {p°,r }, 
where the coordinates of the stochastic vectors p°, r° are determined by formulae 




p(o) = p (o) + jP (o) (G _ 6jP (o) 
p(o) = p (o) + p (o) M + ep (o) 




and 



R { ^=R^+Rf\-d + eR^ 






R 



where ^ = P<°> + P 2 (0) , ^ = R? } + R$\ 
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The next step exactly corresponds to the conflict interaction between systems. 
We introduce new stochastic vectors {p 1 ,!* 1 } with coordinates: 

_(i) = Pj ar j > _(i) _ r j [t a Pj I ■ _ , 9 

^ 1 V-2 (0) (0)''i n (0) (o)' J 

Finally, we have to come back to the non-normalized vectors, which character- 
ize quantitatively populations in both regions after inner and outer conflicts oper- 
ations. So, at time n = 1 we have the following vectors jV _1 {p 1 , r 1 } = {P 1 , R 1 }, 
where 

P 1 = (P 1 (1) ,P 2 (1) ),R 1 = (M 1) ,4 1) ), 



and where 



p(l) _ Jl)^) _ r (l) ? (0) ■ _ i 9 

3 ~ "j P ' 3 ~ 3 R I -* ~ ' 



We can repeat this procedure starting from {P^R 1 }. So we get {P 2 ,R 2 }. 
And so on for any nth step. 

To find the equilibrium points in the case of the complex conflict interac- 
tion described above, we have to solve the following system of equations for 

-Pi) Pi, Ri, Ri- 

(a + l-bR!- cP 1 )(Z 2 - aR 2 (-d + 1 + eP 2 - fR 2 ))Z l = Z, 
\-d+l + ePi - fRi)(Z 2 - aP 2 (a + l-bR 2 - cP 2 ))Z 1 = Z, 
(a + l-bR 2 - cP 2 )(Z 1 - aRii-d + 1 + eP x - fRi))Z 2 = Z, 
\-d+l + eP 2 - fR 2 )(Z 1 -aPi(o + l- W?i - cP 1 ))Z 2 = Z, 

where 

Z x = P x ( a + l-bR 1 - cPi) + Pi(-d + 1 + eP x - fR x ), 
Z 2 = P 2 (a + 1 - bR 2 - cP 2 ) + R 2 (-d + 1 + eP 2 - /P 2 ), 
Z = Z X Z 2 - a[PiP 2 (a + l-bR 2 - cP 2 )(a + 1 - bR x - cP 1 ) + 
+R 1 R 2 {-d + i + e P 1 - fRi){-d + 1 + eP 2 - /P 2 )]. 

We note that in the case a = we have two copies of pure Lotka-Volterra 
models and the corresponding system of equations has at least three equilibrium 
points (trivial, axial, inner positive). 

For the case a ^ it is difficult to obtain exact solutions. Let us obtain some 
insights by numerical approximation. 

Partially, we found that there exist equilibrium points and the limit cycles for 
a wide set of parameter values and initial data (see Figure 4-7). 

Moreover, we established the shift effect for the equilibrium point. Namely, 
we observe that the inner positive equilibrium point (it exist in any system and 
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Figure 6: The existence of the stable oscillations of the conflict interaction 
between Lotka-Volterra systems after 70000 steps of iteration. 
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may be found by formula (4.1)) is shifted after the application of the conflict 
interaction between systems. We see by (4.1) that stabilization of discrete Lotka- 
Volterra model with parameters a = 0.2, b = 0.006, c = 0.002, d = 0.008, e = 
0.002, / = 0, a = 0.007, pf 0) = 3, P 2 (0) = 10, Pj 0) = 5, = 20 occurs when 
Pi = 4, P 2 = 32. This may be easily verified by putting these initial data into 
corresponding equations. In this case we have trivial dynamics. 

Let us consider the case of discrete Lotka-Volterra model with the conflict 
interaction between systems. We take the values of the coefficients a = 0.2, b = 
0.006, c = 0.002, d = 0.008, e = 0.002, / = 0, a = 0.005. Now the equilibrium 
point has the coordinates Pi = 4.043507, P2 = 32.100629. The dynamics is con- 
stant with these initial data. 

In case of larger a, when oscillations appear, the equilibrium point may also 
be easily found if we put the initial data in both systems to be equal. In this case 
the behavior is like in the case of a pure Lotka-Volterra model, and stabilization 
occurs. However, the stable point is shifted, for example, when a = 0.01 (see 
Figure 4-7) the equilibrium point is P 1 = R 1 = 4.087615, P 2 = R 2 = 32.200863. 

Thus, if we have some pray-predator system and want to change the pop- 
ulation inside this system, we may create an analogous "artificial" system, in- 
troduce the conflict interaction and obtain the desired shift of the equilibrium 
point. Apparently a stronger shift of the stable equilibrium point in the case of 
an "ensemble" of larger amount of Lotka-Volterra systems. So, we observed the 
interesting phenomenon: the equilibrium point of an isolated system is shifted if 
we come to the case when identical systems are united as an "ensemble". 

However, this equilibrium point is unstable, any perturbation of initial data 
causes the receding of the system from the equilibrium point. 

One of more interesting observations concerns the limit cycles. It is known 
that no such kind of orbits in discrete Lotka-Volterra model is possible. But 
under the effect of the outer conflict, as we see at the pictures, the dynamical 
system reaches the limit cycle starting both from an inside or outside point with 
respect to the orbit. Partially, in Figure 10,11 we present the model, that starts 
at = 4, P 2 (0) = 32. As it was pointed above, in case of a pure Lotka-Volterra 
model, with these initial data there is no dynamics. However, in the case of the 
model with the outer conflict the process tends to a limit cycle. 

5 Interpretation 

In many works on mathematical biology and economics [SI El El Ell E21 ESI ESI 
ESI ED] the modelling of population dynamics or economical processes is based on 
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Figure 8: Conflict interaction between Lotka- Volterra systems. The starting 
parameters are the same as in Figure 4- 
Initial data are inside the limit cycle (Figure 7). 
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Figure 9: The corresponding phase-space (P x , P 2 )■ Trajectory tends to the 

limit cycle, which is an attractor. 
Initial data are inside the limit cycle (Figure 7). 
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Lotka-Volterra equations. As a rule, continuous, not discrete, models are studied. 
In some works the migration process is considered. It takes place between different 
regions, inside which an interaction of the Lotka-Volterra type is present. For 
example, in [5] the migration rate between regions has some fixed probability. 

We study discrete Lotka-Volterra models with an additional interaction be- 
tween them. That may be interpreted as a some kind of correlation between the 
habitants of different regions. We suppose that discrete models are more natu- 
ral, partially it is clear that birth and death of individuals happen at some fixed 
moments of time. 

It is well known that in the classical discrete pray-predator model a stable 
point exist. The amount of prays and predators tends to this point in the phase- 
space. In this case we observe the following dynamics, after several period of 
oscillations the populations stabilize (see Figure 2). Thus, we have an attracting 
point in phase-space. Such a dynamics exists inside every region when "migra- 
tion" is absent. 

When we introduce an additional interaction between the habitants of differ- 
ent regions a redistribution process appears which we interpret as a migration. 
In some of our complex models there is no stable point, the amount of prays and 
predators in both regions oscillates along fixed orbits. Appearently these orbits in 
a phase-space are attractors. 

We note that explicit formulas of conflict interaction between non-annihilating 
opponents which describe the redistribution of populations are given by (3.1). 
The individuals of a certain kind migrate to the region, where their amount more 
numerous. 

Is the "migration strategy" which is described in our model a natural one? 
We suppose that in many cases individuals may be right behaving in such a way. 
If we consider a pray-predator model, it is clear that every separated individual is 
unable to estimate all factors that have an influence on the population dynamics 
like vital resources inside region, real amount of own and alternative population, 
current population dynamics. In other words, the individual "does not know" 
the parameters of the Lotka-Volterra equations and their current influence on the 
population dynamics. 

However the individual has the group reflex and will migrate to the region, 
where, as he supposes, the vital conditions are best (his population should be 
concentrated there). He suggests, right there are the resources, possibilities for 
reproduction, better conditions to organize large groups. Formula (3.1) just de- 
scribes this tendency. 

Similar motivations may be proposed in case of the work migration. Here the 
unemployed may be regarded as playing the role of " prays" , employees as playing 
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the role of "predators" . People, who seek for work and migrate to another country, 
do not know, as a rule, the real situation in the opposite region. They prefer to 
migrate to the country where the majority of their friends migrated (group reflex). 

Similar, but opposite picture happens with employees who inverts their capital 
to the region with a higher profit. 

So, at the cost of migration accelerates the increasing of one of species' pop- 
ulation in one of the regions. But at the same time there is an effect of the inner 
pray-predator "fight" inside every system. Partially, the population influences 
itself at the cost of the last term in Lotka-Volterra equation. As a result, some 
time later the backward migration starts. 

In the Figure 6 we may see the effect of delay, when the amount of prays 
inside the region decreases, but the predators continue migration to this region, 
until their amount starts decreasing by following the Lotka-Volterra model. 

We emphasize, that in our model, in comparison with discrete Lotka-Volterra 
model, a cyclic oscillations of populations are observed. Moreover, a cyclic at- 
tractor exists in the phase-space, and the pray-predator trajectory tends to this 
orbit both from inside or outside point with respect to this cycle (Figure 8,9). 

We remark that in our model the normalization was fulfilled by the amount 
of habitants of the region, so the component of the corresponding vector may be 
large both at the cost of large population of fixed individuals and at the cost of 
small whole population of the region. So, a migration to the region with a lot of 
"free space" is also possible. 

We also studied model with the attracting interaction (a < 0). In this case we 
obtained formally a similar dynamics, but now individuals migrate to the region 
where they are less numerous. Such a migration strategy might be also natural 
for some species, e.g. for individuals who hunt separately, control large territory 
and have confrontation with relatives. 
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